Historic,  Archive  Document 

Do  not  assume  content  reflects  current 
scientific  knowledge,  policies,  or  practices. 


ft  77*7 

Determination  of  Potential  Direct 
Beam  Solar  Irradiance 


Merrill  R.  Kaufmann  and  James  D.  Weatherred 


Research  Paper  RM-242 
Rocky  Mountain  Forest  and 
Range  Experiment  Station 
Forest  Service 
U.S.  Department  of  Agriculture 


Abstract 


Procedures  are  presented  for  calculating  potential  direct  beam 
solar  irradiance,  corrected  for  latitude,  azimuth  and  inclination  of 
slope,  date,  and  time  of  day.  Equations  are  structured  to  permit  the 
user  to  calculate  instantaneous  or  total  daily  irradiance  using  total 
incoming  shortwave  irradiance  or  any  selected  portion  of  the  solar 
spectrum. 


Kaufmann,  Merrill  R.,  and  James  D.  Weatherred.  1982.  Determina- 
tion of  potential  direct  beam  solar  irradiance.  USDA  Forest  Serv- 
ice Research  Paper  RM-242,  23  p.  Rocky  Mountain  Forest  and 
Range  Experiment  Station,  Fort  Collins,  Colo. 

Procedures  are  presented  for  calculating  potential  direct  beam 
solar  irradiance,  corrected  for  latitude,  azimuth  and  inclination  of 
slope,  date,  and  time  of  day.  Equations  are  structured  to  permit 
the  user  to  calculate  instantaneous  or  total  daily  irradiance  using 
total  incoming  shortwave  irradiance  or  any  selected  portion  of  the 
solar  spectrum. 

Keywords:  Solar  irradiance,  direct  beam  irradiance,  latitude, 
declination,  azimuth,  inclination 


USDA  Forest  Service 
Research  Paper  RM-242 


October  1982 


Determination  of  Potential  Direct 
Beam  Solar  Irradiance 

Merrill  R.  Kaufmann,  Principal  Plant  Physiologist 
Rocky  Mountain  Forest  and  Range  Experiment  Station1 

and 

James  D.  Weatherred,  Computer  Programmer 
Watershed  Systems  Development2 


'Headquarters  is  in  Fort  Collins,  in  cooperation  with  Colorado  State  University. 

2 Fort  Collins  unit  of  USDA  Forest  Service  National  Forest  System,  with  headquarters  in  Washington,  D.C. 


Contents 


Page 


Introduction   1 

Abbreviations   1 

Equations    2 

Horizontal  Surfaces   2 

Tilted  Surfaces   3 

Double  Sunrise  and  Sunset   4 

Obstruction  of  the  Horizon   4 

Explanation  of  Irradiance  Figures   5 

Literature  Cited   5 

Appendix  A    15 

Appendix  B    17 


Determination  of  Potential  Direct  Beam  Solar  Irradiance 

Merrill  R.  Kaufmann  and  James  D.  Weatherred 


Introduction 

Potential  direct  beam  solar  irradiance  is  the  radia- 
tion received  above  the  earth's  atmosphere.  The  max- 
imum direct  beam  irradiance  is  received  when  a  sur- 
face above  the  atmosphere  is  perpendicular  to  the 
sun's  rays.  For  incoming  shortwave  irradiance,  this  is 
equivalent  to  the  solar  constant,  generally  taken  to  be 
1360  W  •  m-2  (1.95  cal  •  cm'2  •  min"1). 

Actual  irradiance  received  at  the  earth's  surface  is 
lower  than  the  potential  direct  beam  irradiance.  The 
reduction  results  from  a  number  of  factors,  some  of 
which  are  highly  predictable  because  they  are  geo- 
metric in  nature,  and  some  of  which  are  less  predict- 
able because  of  varying  atmospheric  effects. 

Geometric  effects  include  latitude,  solar  declination, 
distance  between  the  sun  and  earth  (radius  vector),  in- 
clination and  aspect  of  the  plane  in  question,  time  of 
day,  and  obstructions.  The  plane  or  surface  may  be  a 
hillside,  leaf,  wall  or  roof  of  a  building,  solar  collector, 
etc.  Atmospheric  effects  include  depth  of  clear  at- 
mosphere through  which  radiation  must  pass  (a  func- 
tion of  solar  declination,  latitude,  elevation,  and  time  of 
day),  radiation  scattering  (a  function  of  atmospheric 
composition  and  contamination),  and  cloudiness  (Gates 
1980). 

While  knowledge  of  the  actual  irradiance  received 
at  a  surface  is  useful  to  hydrologists,  plant  physiolo- 
gists, solar  energy  engineers,  and  others,  the  complexi- 
ties of  predicting  actual  irradiance  are  great  because 
atmospheric  conditions  are  dynamic.  Knowledge  of 
potential  direct  beam  irradiance  is  useful  because 
potential  irradiance  is  a  maximum  limit  for  actual  irra- 
diance and  is  usually  used  to  derive  estimates  of  actual 
irradiance. 

Incoming  shortwave  irradiance  is  of  interest  to 
hydrologists,  micrometeorologists,  and  solar  engineers 
because  of  its  direct  role  in  determining  the  energy 
balance  of  a  given  site.  In  some  cases,  however,  only  a 
portion  of  the  solar  spectrum  is  of  concern.  For  exam- 
ple, a  plant  physiologist  studying  photosynthesis  or 
stomatal  behavior  may  be  concerned  only  with  photo- 
synthetic  photon  flux  density  (visible  irradiance,  400  to 
700  nm),  while  a  human  pathologist  evaluating  skin 
disorders  may  be  interested  in  the  ultraviolet  portion  of 
the  spectrum. 

Frank  and  Lee  (1966)  provided  tables  giving  the 
times  of  sunrise  and  sunset  and  the  total  daily  incoming 
irradiance,  which  covered  a  range  of  latitudes, 
aspects,  slopes,  and  times  of  year.  Buffo  et  al.  (1972) 
also  provided  extensive  tables  and  figures  of  direct 
beam  solar  irradiance.  Swift  (1976)  presented  an 
algorithm  for  calculating  daily  total  solar  irradiance  on 
mountain  slopes.  This  algorithm  is  suitable  for  use  as  a 
computer  subroutine,  making  the  use  of  tables  unneces- 


sary. Formulas  for  calculating  direct  beam  irradiance 
presented  by  Frank  and  Lee  (1966),  Gates  (1980),  and 
Lee  (1978)  are  useful,  but  they  require  more  extensive 
development  by  the  user  for  routine  prediction  of  irra- 
diance. Furthermore,  the  formulas  are  considered 
mainly  in  the  context  of  total  incoming  shortwave  irra- 
diance. 

This  paper  has  two  objectives.  First,  the  formulas 
used  by  Frank  and  Lee  (1966),  Buffo  et  al.  (1972),  and 
Swift  (1976)  are  extended  to  provide  more  details  of  the 
calculation  procedures,  particularly  with  regard  to 
determinations  of  hour  angles  for  calculating  sunrise, 
sunset,  and  radiation  and  of  the  effects  of  intervening 
obstruction  by  terrain.  This  is  done  to  facilitate  the  use 
of  computers  for  calculations  of  both  instantaneous 
and  total  direct  beam  irradiance.  The  reader  is  re- 
ferred to  the  programs  listed  in  the  appendixes  and  to 
Swift's  (1976)  algorithm  for  linking  the  calculations.3 

Second,  the  formulas  are  presented  as  a  general 
case  for  direct  beam  solar  irradiance,  where  the  in- 
coming irradiance  may  be  total  shortwave,  visible, 
ultraviolet,  etc.,  as  selected  by  the  user.  In  effect,  the 
formulas  for  instantaneous  and  total  irradiance  yield  a 
"multiplication  factor"  to  be  taken  times  the  potential 
irradiance  (e.g.,  the  solar  constant  for  total  solar  short- 
wave irradiance,  approximately  2,600  fxE  •  m-2  • 
sec-1  for  400  to  700  nm  irradiance,  etc.).  Duffie  and 
Beckman  (1974)  provide  data  from  several  sources 
describing  the  standard  spectral  distribution  of  ex- 
traterrestrial irradiance.  It  must  be  recognized  that 
the  calculation  procedures  given  below  pertain  only  to 
the  determination  of  irradiance  above  the  atmosphere. 
Atmospheric  effects  on  actual  irradiance  received  at  a 
surface  in  question  must  be  determined  for  the  ap- 
propriate portion  of  the  solar  spectrum  of  interest. 

Abbreviations 

MQ       potential  irradiance  constant  (given  a  dimen- 
sionless  value  of  1.0  for  a  surface  above  the 
atmosphere,    perpendicular    to    the  sun's 
rays,  at  the  equinox) 
R0        potential  irradiance  constant  in  user's  units 
Mt        multiplier  for  instantaneous  irradiance  for  a 
given  latitude,  aspect,  inclination,  time  of 
year,  and  time  of  day 
Rj         instantaneous  irradiance  in  user's  units 
Mt        multiplier    for   total    daily   irradiance  from 
sunrise   to    sunset    for    a    given  latitude, 
aspect,  inclination,  and  time  of  year 

3 For  those  who  have  access  to  a  Hewlett  Packard  41 C  hand 
calculator,  a  Users'  Library  Solutions  Manual  for  Solar  Engineer- 
ing (No.  00041-90138)  has  a  program  for  solar-beam  irradiation. 
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Rt         total  daily  irradiance  in  user's  units 

r  radius  vector,  the  ratio  of  the  distance  be- 

tween the  sun  and  earth  at  a  given  date  to 
the  mean  distance 

n         Julian  date 

m         month  of  year  (1  to  12) 

d         day  of  month 

6  latitude  in  degrees  of  the  given  slope,  positive 

in  northern  hemisphere,  negative  in 
southern  hemisphere 

0 '         latitude  of  the  equivalent  slope 

8  solar  declination  in  degrees  for  the  given  slope 

a  change  in  hour  angle  in  degrees  from  the  given 

to  the  equivalent  slope 

co         angular    velocity    of   the    earth's  rotation 
(15°  •  hr"1) 

t  time  in  hours  from  solar  noon 

ta    sunrise  at  the  given  slope 
t2    sunset  at  the  given  slope 
t-[    sunrise  at  the  equivalent  slope 
\'2    sunset  at  the  equivalent  slope 
t"    time  at   the   equivalent   slope  for 
calculating  instantaneous  irradiance 

a  azimuth  of  the  slope  in  degrees  measured 

clockwise  from  north 

i  inclination  of  the  slope,  0°  to  90° 

A         solar  altitude  in  degrees 

A'        solar  azimuth  in  degrees  measured  clockwise 
from  south 


Equations 

Horizontal  Surfaces 

The  basic  equation  for  calculating  instantaneous  ir- 
radiance on  a  horizontal  surface  outside  the  atmos- 
phere is  given  by  Frank  and  Lee  (1966): 

RD 

=  —  (sin0  •  sin  5  +  cos  6  •  cos  5  •  cos  cot).  [1] 
r2 

The  value  of  r2  ranges  from  0.96676  on  January  3  to 
1.03370  on  July  5  (List  1971,  table  169).  As  an  approx- 
imation, r2  can  be  estimated  by  first  determining  the 
Julian  date  and  declination  as  follows: 


By  these  procedures,  8  is  estimated  within  +0.73°  to 
-1.38°,  and  r2  is  estimated  within  +0.84%  to 
-0.70%,  depending  upon  time  of  year. 

Multipliers  for  instantaneous  irradiance  may  be 
calculated  with  a  formula  similar  to  equation  1: 


Mj  =  — —  (sin  6  ■  sin  8  +  cos  6  •  cos  <5  •  cos  cot)  [5] 
r2 

where  MQ  is  assigned  a  value  of  1.0  at  the  equinox.  In 
later  illustrations,  values  of  Mj  are  plotted  as  a  func- 
tion of  several  variables.  The  user  may  calculate  in- 
stantaneous irradiance  using  a  selected  R0  as  follows: 

Rs  =  M,  RQ.  [6] 

For  a  horizontal  surface,  if  the  sum  of  the  absolute 
values  of  6  and  8  is  less  than  or  equal  to  90°,  the  hour 
angles  of  sunrise  (  —  cot)  and  sunset  (cot)  before  and  after 
solar  noon  can  be  found  as  follows: 


cot  =  cos_1(-tan0  ■  tan 5).  [7] 

However,  if  the  sum  of  the  absolute  values  of  6  and  8  is 
greater  than  90°,  then  the  sun  remains  either  above  or 
below  the  horizon  all  day.  If  6  and  8  have  the  same  sign, 
the  sun  remains  above  the  horizon  and  cot  equals  180°. 
If  6  and  8  have  different  signs,  the  sun  is  below  the 
horizon  for  the  day  and  cot  =  0°. 

Times  for  sunrise  (ta)  and  sunset  (t2)  in  hours  before 
or  after  mean  true  solar  noon  are  given  by: 

t2  =  -  tt  =  cot/15.  [8] 

The  solar  times  calculated  above  are  mean  solar 
times.  True  solar  time  differs  from  mean  solar  time  by 
an  amount  that  varies  through  the  year.  True  solar 
time  may  be  calculated  by  algebraically  adding  a  cor- 
rection, called  the  equation  of  time,  which  ranges  from 
- 14  to  + 16  minutes  (List  1971,  table  169).  Frank  and 
Lee  (1966)  give  procedures  for  converting  solar  time  to 
local  standard  time. 

Multipliers  for  total  irradiance  for  a  solar  day  may 
be  calculated  by  integrating  equation  [5]: 


n  =  31(m  -  1)  +  d  -  0.4  m  -  1.8.  [2] 

Round  n  to  the  nearest  day.  If  m  is  1  (January),  add  2  to 
n.  If  m  is  2  (February),  add  3  to  n.  In  leap  years,  add  1 
day  if  m  is  3  or  greater  (Ball  1978). 
Solar  declination  can  be  estimated: 

8  =  23.5  •  sin  [0.9863(284  +  n)].  [3] 

Finally,  r2  is  calculated: 

r2  =  0.999847  +  0.001406(5).  [4] 


Mt  =  (M0/r2)  [2t  •  sin  0  •  sin  8 

+  (12/ir)  •  cos  0  ■  cos  8  •  2  •  sin  cot].  [9] 

The  user  may  calculate  total  irradiance  using  a 
selected  R0  as  follows: 

Rt  =  Mt  •  R0.  [10] 

Note  that  Mt  has  time  units  of  hours;  therefore  R0 
should  have  units  of  hours. 

An  example  set  of  calculations  for  a  horizontal  sur- 
face is  given  in  the  Appendix  (example  1). 
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Tilted  Surfaces 

When  a  surface  is  not  horizontal,  an  "equivalent" 
horizontal  surface  exists  at  a  higher  or  lower  latitude 
and  longitude  which  is  parallel  to  the  given  surface. 
This  equivalent  surface  is  used  to  calculate  the  instan- 
taneous and  total  irradiance  received  on  the  given  sur- 
face. The  next  series  of  formulas,  extending  those  of 
Frank  and  Lee  (1966,  citing  earlier  references),  deter- 
mine the  location  of  the  equivalent  surface.  Refer  to 
figure  1  for  a  depiction  of  a  given  and  equivalent  slope. 

The  latitude  of  the  equivalent  surface  is  determined 
as  follows: 

6'  =  sin-1  (sin  i  •  cos  a  •  cos  0  +  cos  i  •  sin0)  [11] 

where  inclination  is  in  degrees  and  azimuth  is  0°  to 
360°.  If  inclination  is  given  in  percent  slope,  then: 

i  =  tan  - 1  (%  slope/100).  [12] 

The  longitudinal  correction  of  hour  angle  for  the 
equivalent  surface  is  given  by: 

a  =  tan  - 1  [(sin  a  •  sin  i) 

-5-  (cos  i  •  cos  0  -  cos  a  •  sin  i  •  sin  0)]  [13] 

when  a  is  greater  than  0°. 

Because  the  tan-1  of  an  angle  is  defined  for  angles 
between  -90°  and  90°,  if  the  slope  is  east-facing  and  a 
is  less  than  0°,  then  a  is  added  to  180°.  If  the  slope  is 
west-facing  and  a  is  greater  than  0°,  a  is  added 
to  - 180°.  If  a  is  0°  (due  north)  and  the  sum  of  1 0|  and  i  is 


Figure  1.— Relation  of  given  slope  and  equivalent  slope.  Given 
slope  at  latitude  6  has  azimuth  a  and  inclination  i.  Equivalent 
slope  has  a  different  latitude,  6' ,  and  the  longitude  differs  from 
that  of  the  given  slope  by  the  hour  angle,  a. 


less  than  or  equal  to  90°,  a  equals  0°.  However,  if  the 
sum  of  |0|  and  i  is  greater  than  90°,  a  equals  180°. 

Hour  angles  for  sunrise  and  sunset  at  the  equivalent 
surface  are  calculated  the  same  way  as  those  for 
horizontal  slopes  in  equation  [7],  except  that  the 
equivalent  latitude  is  used: 

cot'  =  cos-1(~tan0'  -  tan 5).  [14] 

The  same  limits  apply  to  (-tan  0'  •  tan  5)  and  cot'  using 
6'  as  to  (-tan  0  •  tan  <5)  and  cot  for  equation  [7]. 
The  special  case  of  a  double  sunrise  and  sunset  is  con- 
sidered later. 

Hour  angles  for  sunrise  and  sunset  of  the  given  slope 
may  be  determined  as  follows.  If  (  —  cot  -I-  a)  is  greater 
than  or  equal  to  cot'  (values  from  eq.  [7],  [13],  and  [14]), 
then  the  given  slope  is  completely  shaded  throughout 
the  day  (ta  and  t2  equal  0).  If  the  slope  is  not  completely 
shaded,  the  time  of  sunrise  at  the  given  slope  then 
becomes: 

tt  =  [maximum  of  ( -  cot)  and  ( -  cot '  -  a]]  1 15.  [15] 

The  time  of  sunset  at  the  given  slope  becomes: 

t2  =  [minimum  of  (cot)  and  (cot'  -  a)]  / 15.  [16] 

Irradiance  calculations  for  the  equivalent  surface 
require  sunrise  and  sunset  times  for  the  equivalent  sur- 
face. Sunrise  for  the  equivalent  slope  is: 

t^  =  [maximum  of  ( -  cot  +  a)  and  ( -  cot ' )]  / 15.  [17] 

Sunset  for  the  equivalent  slope  is: 

t2  =  [minimum  of  (cot  +  a)  and  (cot')]  / 15.  [18] 

Finally,  irradiance  multipliers  can  be  calculated.  In- 
stantaneous irradiance  is  determined  as  follows: 

M 

Mj  =  — -  [sin0'  •  sin  5  +  cos0'  •  cos  5  •  cos  (cot")]  [19] 
r2 

where  cot "  =  cot  +  a  and  t  is  greater  than  or  equal  to  ta 
and  less  than  or  equal  to  t2.  If  t  is  less  than  tr  or 
greater  than  t2,  then  Mj  equals  0. 

A  series  of  graphs  discussed  below  present  Mj  as  a 
function  of  latitude,  azimuth,  inclination,  date,  and 
time  of  day.  Recalling  equation  [6],  instantaneous  irra- 
diance in  user  units  (Rj)  can  be  determined  as  the  prod- 
uct of  Mj  and  R0.  Alternatively,  the  numerator  value  of 
M0  in  equation  [19]  may  be  replaced  with  RQ. 

Multipliers  for  total  daily  irradiance  are  calculated: 

Mt  =  (M0/r2)  [(t;-  y  •  sin  0'  •  sin  8 

+  (12/7T)  •  cos  0'  •    COS  8 

•  (sin  cot2  -  sin  cot^)].  [20] 

Total  irradiance  in  user  units  (Rt)  is  determined  as  the 
product  of  Mt  and  RQ,  or  Rt  may  be  determined  by 
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substituting  RQ  for  the  numerator  value  of  M0  in  equa- 
tion [20].  Values  of  Rt,  using  the  solar  constant  as  R0, 
for  a  range  of  latitudes,  azimuths,  inclinations,  and 
dates  are  tabulated  in  Frank  and  Lee  (1966). 

An  example  set  of  computations  for  a  tilted  surface 
is  given  in  the  Appendix  (example  2). 

Double  Sunrise  and  Sunset 

Steep,  north-facing  slopes  and  the  north  sides  of 
buildings  at  middle  latitudes  in  the  northern  hemi- 
sphere exhibit  a  sunrise  and  sunset  both  in  the  morning 
and  in  the  afternoon.  Double  sunrises  and  sunsets  oc- 
cur when  the  following  conditions  are  met: 

(-cot')<  (cot  +  a  -  360) <  (-cot  +  a)  <  (cot').  [21] 

If  these  conditions  exist,  sunrise  and  sunset  times  for 
the  given  slope  are  calculated  as  follows: 


First  tt  =  -  cot  / 15  [22] 

First  t2  =  (cot'  -a)  1 15  [23] 

Second  t,  =  (360  -  cot '  -  a)  1 15  [24] 

Second  t2  =  cot  / 15.  [25] 

For  calculating  irradiance  at  the  equivalent  surface, 
sunrise  and  sunset  times  are: 

First  t[  =  -cot'  /  15  [26] 

First  t;  =  (cot  +  a  -  360)  /  15  [27] 

Second  t[  =  (-cot  +  a)  1 15  [28] 

Second  V2  =  cot'  / 15.  [29] 


Calculation  of  the  total  daily  irradianco  multiplier 
for  days  having  two  sunrises  and  sunsets  requires 
determining  Mt  using  equation  [20]  for  the  first  solar 
day  (first  t{  to  t^)  and  for  the  second  solar  day  (second 
t[  to  t^).  These  values  are  added  to  obtain  the  total 
daily  irradiance  multiplier.  Again,  total  daily  irra- 
diance in  user's  units  may  be  determined  as  the  pro- 
duct of  Mt  and  R0  or  by  replacing  the  numerator  value 
of  MQ  with  RQ  in  equation  [20]. 

An  example  set  of  calculations  for  double  sunrise 
and  double  sunset  times  is  given  in  the  Appendix  (ex- 
ample 3).  For  radiation  calculations  on  days  of  double 
sunrise  and  sunset,  refer  to  example  2. 

Obstruction  of  the  Horizon 

In  mountainous  regions,  few  sites  are  characterized 
by  a  direct  view  of  the  earth's  horizon.  Unless  sites  are 
located  upon  the  highest  ridgetops,  intervening  hills 
and  ridges  obscure  the  sun,  resulting  in  later  sunrises 
or  earlier  sunsets  than  those  predicted  by  equation  [8]. 
The  delay  of  sunrise,  advancement  of  sunset,  or  midday 
obstruction  by  intervening  terrain  (or  objects)  can  be 
determined  by  comparing  the  elevation  angle  of  the 
obstruction  with  the  solar  altitude.  When  the  solar 
altitude  at  a  given  azimuth  is  below  the  obstruction 
altitude  at  the  same  azimuth,  irradiance  is  set  to  zero. 


Procedures  for  determining  the  solar  altitude  and 
azimuth  are  given  by  Lee  (1978).  The  solar  altitude  is 
calculated  as  follows: 

A  =  sin-1  (sin  0  •  sin  5  +  cos  6  •  cos  5  ■  cos  cot).  [30] 

Values  of  A  less  than  0°  indicate  that  the  sun  is  below 
the  observer's  true  horizon,  while  values  greater  than 
0°  indicate  that  the  sun  is  above  the  observer's  true 
horizon. 

The  computation  of  solar  azimuth  is  somewhat  more 
difficult.  Solar  azimuths  calculated  below  are  deter- 
mined clockwise  from  the  observer's  south.  If  |0|  is 
greater  than       the  solar  azimuth  is  given  by 

A'  =  sin-1  (cos  6-  sin  cot  /  cos  A).  [31] 

To  establish  the  correct  quadrant,  first  calculate  the 
hour  angle  (Hw)  at  which  the  sun  is  due  west  of  the 
observer: 

Hw  =  cos  - 1  (tan  5  •  cotan  0).  [32] 

Then,  if  |cot|  is  less  than  or  equal  to  Hw,  A'  is  in  the 
correct  quadrant.  However,  if  |cot|  is  greater  than  Hw, 
A '  is  placed  in  the  correct  quadrant  by 

A'  =  sign (180,  cot)  -  A'  [33] 

where  sign  (180,  cot)  is  180°  if  cot  is  greater  than  or 
equal  to  0  and  -  180°  if  cot  is  less  than  0°. 

If  |0|  is  less  than  or  equal  to  three  possible 
cases  exist.  In  the  first  case,  |cot|  is  greater  than  0° 
but  less  than  180°.  If  8  equals  0°,  then  0  equals  0°  and 

A'  =  sign (90,  cot).  [34] 

If  8  is  not  equal  to  0°,  then  A'  is  calculated  as  in  equa- 
tion [31].  A'  from  equation  [31]  is  in  the  correct 
quadrant  if  8  is  less  than  0°.  For  8  greater  than  0°,  A' 
is  determined  using  equation  [33]. 

In  the  second  case,  cot  equals  0°.  If  |0|  equals 
then  when  0  is  greater  than  or  equal  to  0°,  A'  is  0°; 
when  0  is  less  than  0°,  A'  is  180°.  If  |0|  is  less 
than       then  when  8  is  greater  than  0°,  A'  is  180°; 
when  8  is  less  than  0°,  A'  is  0°. 

In  the  third  case,  |cot|  equals  180°.  If  |0|  equals  |6|, 
then  when  0  is  greater  than  or  equal  to  0°,  A'  is  180°, 
and  when  0  is  less  than  0°,  A'  is  0°.  If  |0|  is  less  than  |5|, 
then  when  8  is  greater  than  0°,  A'  is  180°,  and  when  8 
is  less  than  0°,  A'  is  0°. 

To  determine  the  effects  of  intervening  terrain,  the 
terrain  altitude  and  azimuth  must  be  determined  for 
the  opposite  horizon  (east  horizon  for  a  slope  having  an 
azimuth  of  0°  to  180°  and  west  horizon  for  a  slope  hav- 
ing an  azimuth  of  180°  to  360°).  This  may  be  done  by 
mapping  the  elevation  angle  of  the  horizon  at  suitably 
spaced  azimuths  or  by  calculating  the  elevation  angle 
using  topography  maps. 
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At  some  sites,  such  as  lower,  south-facing  slopes  in 
valleys  (northern  hemisphere),  the  entire  horizon  be- 
tween the  points  of  sunrise  and  sunset  must  be  mapped 
to  determine  if  the  sun  is  obstructed  by  the  opposite 
ridge.  The  portion  of  the  obstructed  horizon  which 
must  be  mapped  and  the  azimuth  difference  between 
points  at  which  elevation  angle  is  measured  will  de- 
pend upon  the  particular  site  in  question  and  the  ac- 
curacy required. 

After  the  terrain  profile  is  determined,  solar  altitude 
is  compared  with  the  terrain  elevation  at  various  solar 
azimuths  to  determine  the  times  at  which  the  obstruc- 
tion occludes  the  sun.  For  instantaneous  irradiance 
calculations,  Mi  or  Rj  is  set  to  zero  if  the  solar  altitude 
is  less  than  the  elevation  angle  of  the  terrain.  For  total 
daily  irradiance,  Mt  or  Rt  are  calculated  using  t[  and  \.'2 
adjusted  for  the  effects  of  slope  at  the  given  site  and  for 
the  effects  of  terrain  elevation. 


Explanation  of  Irradiance  Figures 

A  series  of  graphs  was  prepared  to  illustrate  how 
the  instantaneous  irradiance  multiplier  is  influenced 
by  latitude,  date,  azimuth  (aspect),  inclination  (percent 
slope),  and  time  of  day.  Three  figures  are  presented  for 
each  of  three  latitudes— 30°,  40°,  and  50°  N.  These 
figures  are  for  June  22,  March  or  September  22,  and 
December  22,  representing  the  summer  solstice,  vernal 
and  autumnal  equinoxes,  and  winter  solstice  in  the 
northern  hemisphere.  Slopes  in  degrees  are  as  follows: 
0%  (0°),  20%  (11.3°),  50%  (26.6°),  100%  (45°),  200% 
(63.4°),  400%  (76.0°).  Effects  of  obstructed  opposite 
horizons  are  not  taken  into  account. 

The  figures  illustrate  the  complex  geometric  rela- 
tionships among  latitude,  solar  declination,  azimuth,  in- 
clination, and  hour  angle.  Total  day  length  decreases 
from  June  22  to  December  22,  but  the  change  in  day 
length  is  greatest  at  high  latitudes.  For  example,  at  30° 
N  latitude,  day  length  decreases  from  13.9  hours  to 
10.1  hours  (figs.  2  and  4),  while  at  50°  N  latitude,  day 
length  decreases  from  16.2  hours  to  7.8  hours  (figs.  8 
and  10). 

Vertical  dashed  lines  in  the  figures  represent 
sunrises  or  sunsets  at  the  true  horizon  for  surfaces 
tilted  east  or  west.  For  east-facing  surfaces  (azimuth  of 
90°),  an  increase  in  slope  increases  the  intensity  of 
radiation  immediately  after  the  sun  rises  past  the 
horizon.  However,  steeper  slopes  also  experience 
earlier  sunsets  because  the  effective  horizon  is 
elevated.  Opposite  effects  occur  on  west-facing  slopes. 

Instantaneous  irradiance  of  north-facing  slopes 
(azimuth  of  0°)  at  solar  noon  decreases  as  the  inclina- 
tion increases.  At  steep  inclinations,  the  sun  may  set  in 
the  morning  and  rise  in  the  afternoon  during  summer 
months;  this  effect  becomes  more  pronounced  at  higher 
latitudes  (compare  figs.  2,  5,  8).  During  winter  months, 
low  sun  angles  prevent  any  direct  beam  radiation  from 
reaching  steep,  north-facing  slopes  (figs.  4,  7,  10).  On 
south-facing  slopes  (azimuth  of  180°),  instantaneous  ir- 
radiance is  often  higher  than  on  horizontal  surfaces, 


because  the  normal  to  the  slope  is  close  to  the  solar 
beam  pathway. 

Figures  2  through  10  are  presented  only  to  demon- 
strate the  effects  of  various  factors  which  influence  in- 
stantaneous irradiance.  Through  proper  use  of  the 
equations  given  above,  the  user  may  calculate  instan- 
taneous irradiance  for  any  latitude,  date,  azimuth,  or 
inclination.  In  the  southern  hemisphere,  instantaneous 
irradiance  multipliers  for  slopes  facing  directly  east  or 
west  are  similar  to  those  shown  in  the  figures. 
However,  for  north-  and  south-facing  slopes,  the 
graphs  are  reversed. 

The  total  irradiance  multiplier  (Mt),  calculated  with 
appropriate  equations  above,  is  represented  by  the  area 
under  the  curves  between  sunrise  and  sunset.  Frank  and 
Lee  (1966)  tabulated  total  daily  irradiance  for  latitudes 
between  30°  and  50°  N,  16  azimuths,  and  slopes  from  0% 
to  100%  for  various  times  of  the  year.  Their  values  were 
calculated  with  a  solar  irradiance  constant  of 
2.0  cal  •  cm-2  •  min-1.  Irradiance  multipliers  for  total 
solar  irradiance  (Mt)  are  equal  to  Frank  and  Lee's  values 
divided  by  120  (from  2  cal  •  cm-2  •  min-1  times  60  min). 
Tabulated  values  of  Buffo  et  al.  (1972)  include  an  at- 
mospheric transmission  coefficient  and  cannot  be  com- 
pared directly  with  those  given  here  or  in  Frank  and  Lee 
(1966). 
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Figure  2.— Instantaneous  irradiance  multiplier  for  30°  N  latitude  on  June  22. 
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Figure  3.— Instantaneous  irradiance  multiplier  for  30°  N  latitude  on  March  22  or  September  22. 
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Figure  4.— Instantaneous  irradiance  multiplier  for  30°  N  latitude  on  December  22. 
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Figure  5.— Instantaneous  irradiance  multiplier  for  40°  N  latitude  on  June  22. 
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Figure  6.— Instantaneous  irradiance  multiplier  for  40°  N  latitude  on  March  22  and  September  22. 
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Figure  7.— Instantaneous  irradiance  multiplier  for  40°  N  latitude  on  December  22. 
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Figure  8.— Instantaneous  irradiance  multiplier  for  50°  N  latitude  on  June  22. 
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Figure  9.— Instantaneous  irradiance  multiplier  for  50°  N  latitude  on  March  22  or  September  22. 


13 


Latitude:     50  deg  N 


Date:       December  22 


1.2  + 


Azimuth:  0  deg  %  Slope 


0.8 


0.4  ■• 


u 

a 


22  0.0 
c     1.2  ♦ 


4J 


o  o 
o  o 
o  o 

o  o 
o  o 
ox        X  o 

X  X 


T3 
03 
05 


0.8  - 


Azimuth:  180  deg 

•Mi 
§  1 
I  +++  I 

f  *  ♦! 

#  » 
l#     .*  ? 


+  • 


.  + 


0.4  ■• 


x*x 

X  X 
X  X 
X  X 


X      o      O  x 
o  o 
o  o 

'X  X 

o  o 


0.0 


12 


o 

X 

# 


0 
20 
50 
100 
200 
400 


Azimuth:   90  deg 


...t 


x  o'T^o 
X    o  $  *  o 

x    °    a+*X  ° 
O  X  o 

O  *+•  X 

o  1     *  X 

■  #+ 


Azimuth:   270  deg 


i 

qJ?*o  x 

o  *•  +  o 
ox  o 

O  X  +r 


o  X  ' 


+4 


24  0 


12 


Time  of  Day  (hr) 


Figure  10.— Instantaneous  irradiance  multiplier  for  50°  N  latitude  on  December  22. 


14 


Appendix  A 


Example  1. — Horizontal  Surface 

Date  June  22 

6  40° 

5  23.5° 

R0  1360  W  •  m"2  (1.95  cal  ■  cm"2   •  min"1) 

t  -2.00  hr  (1000  hr) 

Times  of  Sunrise  and  Sunset: 

cot    ^cos"1  (-tan  40°  •  tan  23.5°) 
t2    =  7.43  hr  (1926  hr) 
tt    =  -7.43  hr  (0434  hr) 

Instantaneous  Irradiance  Multiplier: 

cot    =  -2  •  15°  =  -30° 
Mj  =1.0[sin40°  ■  sin  23.5° 

+  cos  40°  •  cos  23.5°  ■  cos  (-30°)] 

-  (0.999847  +  0.001406  ■  23.5°) 

=0.8372 

Instantaneous  Irradiance: 

R;    =0.8372  •  1360  W  •  m"2 
Rj    =1138.6  W  ■  m~2 

Total  Daily  Irradiance  Multiplier: 
Mj  =1.0[2  •  7.43  •  sin  40°  ■  sin  23.5° 

+  (12/tt)  •  cos  40°  •  cos  23.5°  •  2  •  sin  111.4°] 

-  1.0329 
M,  =8.5251 

Total  Daily  Irradiance: 

Rt    =8.5251  •  60  •  1360 
Rt    =6.956  X  105  W  •  m"2 


Example  2. — Tilted  Surface 


Date 

September  22 

e 

40° 

8 

0° 

1360  W  •  m"2 

t 

1.00  hr  (1300  hr) 

a 

90° 

i 

40% 

Inclination  in  Degrees: 

i  =  tan"1  (40/100) 
i  =  21.8° 

Equivalent  Latitude: 

6'  =  sin'1  (sin  21.8°  •  cos  90°  •  cos  40° 
+  cos  21.8°  •  sin  40°) 

6'  =  36.64° 


Hour  Angle  Correction: 

a  =  tan-1  [(sin  90°  •  sin  21.8°) 

-s-  (cos  21.8°  •  cos  40° 

-  cos  90°  •  sin  21.8°  ■  sin  40°)] 

a  =  27.57° 


Sunrise  and  Sunset  Hour  Angles  on  Equivalent  Slope: 

wt'=  cos"1  (-tan  36.64°  •  tan  0°) 

cot'  =  90° 

Also, 

cot  =  cos"1  (-tan  40°  •  tan  0°) 
cot  =  90° 

Sunrise  and  Sunset  for  the  Given  Slope: 

t1  =  [max(-90°,  -  117.57°)]  / 15 

t1  =  -  6.00  hr  (0600  hr) 

t2  =  [min(90°,  62.43°)]  / 15 

t2  =  4.16  hr  (1610  hr) 

Sunrise  and  Sunset  for  the  Equivalent  Slope: 

t{  =  [max(- 62.43°,  -90°)]/ 15 
tj  =  -4.16  hr  (0750  hr) 
^  =  [min(117.57°,  90°)]  / 15 
V2  =  6.00  hr  (1800  hr) 

Instantaneous  Irradiance  Multiplier: 

cot"  =  15  •  1.00  +  27.57° 

cot"  =  42.57° 

Mi  =  1.0[(sin  36.64°  •  sin0°) 

+  (cos  36.64°  •  cos  0°  ■  cos  42.57°)] 
-v-  (0.999847  +  0.001406  ■  0°) 
Mj  =  0.5910 

Instantaneous  Irradiance: 

Rj  =  0.5910  •  1360  W  ■  m~2 
Rj  =  803.8  W  •  m"2 

Total  Daily  Irradiance  Multiplier: 

Mt  =1.0{(6.00  +  4.16)  •  sin  36.64°  ■  sin  0° 
+  (12/tt)  •  cos  36.64°  ■  cos  0° 
•  [sin  90°   -  sin  (-62.43°)]} 
-  (0.999847  +  0.001406  ■  0°) 

Mt  =  5.7827 

Total  Daily  Irradiance: 

Rt  =  5.7827  •  60  ■  1360  W  •  m~2 
Rt  =  4.719  X  105  W  •  m  2 
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Example  3. — Double  Sunrise  and  Sunset 


Date 

June  22 

e 

50° 

d 

23.5° 

a 

0° 

i 

250% 

Inclination  in  Degrees: 

i  =  tan"1  (250/100) 
i  =  68.2° 

Equivalent  Latitude: 
d'  =  sin"1  (sin  68.2°  •  cos  0°  •  cos  50° 

+  cos  68.2°  •  sin  50°) 
6'  =  61.8°  (note:  this  latitude  is  on  the  opposite  side 

of  the  North  Pole  from  the  given  slope;  see 

calculation  for  a  below) 

Hour  Angle  Correction: 

6  +  i  =  50  +  68.2 
6  +  i  =  118.2 

Because  [0  +  i)  >  90°,  a  =  180° 
cot  =  cos"1  (-tan  50°  •  tan  23.5°) 
cot  =  121.21° 


cot'  =  cos-1  (-tan  61.8°  •  tan  23.5°) 
cot'  =  144.19° 

Times  of  Double  Sunrise  and  Sunset  at  the  Given  Slope: 

first  ta  =  -  121.21°/15 

first  t1=  -  8.08  hr  (0355  hr) 

first  t2  =  (144.19°  -  180°)  / 15 

first  t2  =  -  2.39  hr  (0937  hr) 

second  tt  =  (360°  -  144.19°  -  180°)  / 15 

second  ta  =  2.39  hr  (1423  hr) 

second  t2  =  121.21°/15 

second  t2  =  8.08  hr  (2005  hr) 

Times  of  Double  Sunrise  and  Sunset  at  the  Equivalent 
Slope: 

first  t[=  - 144.19°/15 

first  t[=  -  9.61  hr  (0223  hr) 

first  t2  =  (121.21°  +  180°  -  360°)  / 15 

first  t2  =  -  3.92  hr  (0805  hr) 

second^  =  (-121.21°  +  180°)  / 15 
second  t[  =  3.92  hr  (1555  hr) 
second  t2  =  144.19°/15 
second  t2  =  9.61  hr  (2137  hr) 
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U.S.  Department  of  Agriculture 
Forest  Service 

Rocky  Mountain  Forest  and 
Range  Experiment  Station 


The  Rocky  Mountain  Station  is  one  of  eight 
regional  experiment  stations,  plus  the  Forest 
Products  Laboratory  and  the  Washington  Office 
Staff,  that  make  up  the  Forest  Service  research 
organization. 

RESEARCH  FOCUS 

Research  programs  at  the  Rocky  Mountain 
Station  are  coordinated  with  area  universities  and 
with  other  institutions.  Many  studies  are 
conducted  on  a  cooperative  basis  to  accelerate 
solutions  to  problems  involving  range,  water, 
wildlife  and  fish  habitat,  human  and  community 
development,  timber,  recreation,  protection,  and 
multiresource  evaluation. 

RESEARCH  LOCATIONS 

Research  Work  Units  of  the  Rocky  Mountain 
Station  are  operated  in  cooperation  with 
universities  in  the  following  cities: 

Albuquerque,  New  Mexico 

Bottineau,  North  Dakota 

Flagstaff,  Arizona 

Fort  Collins,  Colorado* 

Laramie,  Wyoming 

Lincoln,  Nebraska 

Lubbock,  Texas 

Rapid  City,  South  Dakota 

Tempe,  Arizona 


'Station  Headquarters:  240  W.  Prospect  St.,  Fort  Collins,  CO  80526 


